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In this work we have studied the possibility of obtaining cosmic acceleration in Brans-Dicke theory 
with varying or constant uj (Brans- Dicke parameter) and with or without self-interacting potential, 
the background fluid being barotropic fluid or Generalized Chaplygin Gas. Here we take the power 
law form of the scale factor and the scalar field. We show that accelerated expansion can also be 
00 ■ achieved for high values of uj for closed Universe. 
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I. INTRODUCTION 

Recent measurements of redshift and luminosity-distance relations of type la Supernovae indicate that the 
expansion of the Universe is accelerating [1, 2]. This observation gives rise to the search for a matter field, 
which can be responsible for accelerated expansion. There are several proposals regarding this, Cosmological 
Constant, Quintessence, Dark Energy [3 - 5] being some of the competent candidates. However, most of 
these models fit only to spatially fiat (k — 0) Friedmann- Robertson- Walker model [6], though a few models 
[7] work for open Universe (k = —1) also. Brans- Dicke (BD) theory has been proved to be very effective 
_j regarding the recent study of cosmic acceleration [8]. BD theory is explained by a scalar function <f> and a 
constant coupling constant u>, often known as the BD parameter. This can be obtained from general theory 
of relativity (GR) by letting uj — > oo and <f> — constant [9]. This theory has very effectively solved the 
problems of inflation and the early and the late time behaviour of the Universe. N. Banerjee and D. Pavon 
[8] have shown that BD scalar tensor theory can potentially solve the quintessence problem. The generalized 
BD theory [10] is an extension of the original BD theory with a time dependent coupling function u>. In 
Generalized BD theory the BD parameter uj is a function of the scalar field <f>. N. Banerjee and D. Pavon 
have shown that the generalized BD theory can give rise to a decelerating radiation model where the big-bang 
nucleosynthesis scenario is not adversely affected [8] . Modified BD theory with a self-interacting potential have 
also been introduced in this regard. Bertolami and Martins [11] have used this theory to present an accelerated 
Universe for spatially flat model. All these theories conclude that uj should have a low negative value in order 
!* ■ to solve the cosmic acceleration problem. This contradicts the solar system experimental bound uj > 500. 
However Bertolami and Martins [11] have obtained the solution for accelerated expansion with a potential 
4> 2 and large \uj\, although they have not considered the positive energy conditions for the matter and scalar field. 
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In this paper, we investigate the possibilities of obtaining accelerated expansion of the Universe in BD theory 
where we have considered a self-interacting potential V which is a function of the BD scalar field <f> itself and 
a variable BD parameter which is also a function of <fi. We show all the cases of uj — constant, uj = uj(4>), 
V = and V = V((f>) to consider all the possible solutions. We examine these solutions for both barotropic 
fluid and the Generalized Chaplygin Gas [12 - 13], to get a generalized view of the results in the later case. We 
analyze the conditions under which we get a negative q ( deceleration parameter, — in all the models of the 
Universe. For this purpose we have shown the graphical representations of these scenario for further discussion. 

The paper is organized as follows: In section II, the field equations for self-interacting BD theory have been 
given. Sections III and IV deals with the different cases of barotropic fluid and Generalized Chaplygin Gas 
respectively. Each of these two sections are divided into two parts, namely A(V^ = 0) and B (V = V(4>)) where 
again two different cases have been considered with uj = constant and uj = uj(4>) respectively We have taken 
some particular values of the constants for the graphical representations of V and w against the variation of 
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the scalar field <f>. We have discussed the results obtained in section V. 



II. FIELD EQUATIONS 

The self- interacting Brans-Dicke theory is described by the action: (choosing 8irGo = c = 1) 
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V{4>) + Cr, 
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where V((f>) is the self-interacting potential for the BD scalar field <j) and u>(4>) is modified version of the BD 
parameter which is a function of </> [9]. The matter content of the Universe is composed of perfect fluid, 



Tpv = {p + p)u^u v + p g^ u 
where u v = — 1 and p, p are respectively energy density and isotropic pressure. 
From the Lagrangian density (1) we obtain the field equations 
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3 + 2u(<l>) 3 + 2w(0) [ ■ - rr/ T d(j) 
where T = T^g^ . 

The line element for Friedman-Robertson- Walker spacetime is given by 

dr 2 



3 + 2u(<f>y^ 



(3) 



(4) 



ds 2 = -dt A + a 2 (t) 



+ r 2 (d6 2 + sin 2 6d<j> 2 ) 



(5) 



1 - kr 2 

where, a(t) is the scale factor and k(= 0, ±1) is the curvature index. 

The Einstein field equations for the metric (5) and the wave equation for the BD scalar field cf> are the following 



a 2 + k p acf) lj (jr V(4>) 
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The energy conservation equation is 



p + 3-(p + p) = Q 



(7) 



(8) 



(9) 



Now we consider two types of fluids, first one being the barotropic perfect fluid and the second one is 
Generalized Chaplygin gas [12, 13]. 
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III. MODEL USING BAROTROPIC FLUID 

Here we consider the Universe to be filled with barotropic fluid with EOS 

P = 1P (-1<7<1) (10) 

The conservation equation (9) yields the solution for p as, 

p = Po a-^ (11) 

where po(> 0) is an integration constant. 



A. Solution without potential: V{<f>) = 

Case I: 

First we choose oj(4>) = w=constant. 

Now we consider power law form of the scale factor 

a(t) = a t a (a > 1) (12) 
In view of equations (10) and (11), the wave equation leads to the solution for <f) to be 



Poao 



-3(l+7) t 2-3a(l+ 7 ) 



(2w + 3)(l-3o7)[2-3a(l+7)] 
For fc^Owe get from the field equations (6) and (7), the value of a = 1 and 



(13) 



^ (7 -l)(3 7 +!)-!- A 
z do 



(37 + 1) 

We have seen that 7^-^ and we have 

2(1 + A 



( 7 -l)(3 7 + l) 



(14) 



(15) 



Since u> must be negative for — | < 7 < 1, we have seen that for this case the deceleration parameter q = 0, 
i.e., the universe is in a state of uniform expansion. For k = 0, the field equations yield 

[2 - 3a( 7 + 1)] [2(2a - 1) + wfr - 1){2 - 3a( 7 + 1)}] = (16) 
From equation (16) we have two possible solutions for a: 



a 



3(7+1) 



IT for -1 < 7 < -i 



aild a = [4+L"(l-7 7 4 fOT -| < ^ <1 

For these values of a, we have seen that ui < and the deceleration parameter q < 0. Thus for k = 
with the power law form of the scale factor a = a^t a it is possible to get the accelerated expansion of the Universe. 

Case II: 



Now we choose u> = w{4>) to be variable. Here we consider the power law form of <f) as 
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Fig. 1 and 2 shows the variation of to against <f> for different values of a = 1, 1.5, 2, 2.5, 3. In Fig 1 we have considered 
fiat model, i.e., k — and the present dust filled epoch, i.e., 7 = 0, normalizing the parameters as a u = po = 4>a = 1, 
also as the calculation shows, for this f3 = —2, whereas for Fig 2 we have taken closed model, i.e., k = 1 and 
fflo = 0o = 1, po = 6, 7 = i and /3 = — 2a, according to the calculations. 



m = to*" (17) 

with the power law from of a(t) given by equation (12). 
Proceeding as above we get 



ad + 2a + 3 - 3 2 1 + 7 , 3e»( 7 +i)-2 3a ( 7 +i)+/3-2 2fc 2(1-0) 
— al TyTPoao 3(1+7 Vo " 4> » + -TT-T0 " 



Now for acceleration q < implies that a > 1. Using the other equations we arrive at two different situations: 

(i) First considering the flat Universe model, i.e., k = 0, we get, 3=1 — 3a, i.e., 3 < — 2 for 7 > | , 3 = —2a, 
i.e., /3 < — 2 (as a > 1) for 7 = | and 3 = — 2 for 7 < |. That is cosmic acceleration can be explained at all 
the phases of the Universe with different values of 3 where <f) = (f>ot^ 

(ii) If we consider the non-flat model of the Universe, i.e., k ^ 0, we are left with two options. For closed 
model of the Universe, i.e., for k = 1 we can explain cosmic acceleration for the radiation phase only and for 
that 3 = —2a giving 3 < —2 and 6</>o a o 2 = po, whereas we do not get any such possibility for the open model 
of the Universe. 

Now preferably taking into account the recent measurements confirming the flat model of the Universe, if 
3 = — 2 we see that we have an accelerated expansion of the Universe after the radiation period preceded by 
a decelerated expansion before the radiation era and a phase of uniform expansion at the radiation era itself. 
Also if 3 < — 2 cosmic acceleration is followed by a deceleration phase as a < 1 for 7 < |. 



B. Solution with potential: V = V{4>) 

Case I: 

Let us choose u)(<p) = u =constant. 

In this case instead of considering equations (12) and (17) we consider only one power law form 
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<t> - ^a a (19) 
Using equation (19) in equations (6) and (7) we get 



a = 



2k + 2(l + 7 )- P ° 



i> {37a + 67 - a 2 + la + 6 - 2wa 2 } 



Putting fc = 0, we get 

a = At (20) 



p (l+7)(3+a+37) 2 

2po{6(l+7)+a(7+37)-a 2 (l+2w)} 



3+a + 3 7 



where ^4 

Therefore, = Bt 3 +°>"+ 3 ~i where, B = 4> A a 
Now ' if 3+a+3 7 ^ we § et 

a < -(l + 3 7 ) (21) 
Substituting these values in (6), (7), (8), the solution for the potential V is obtained as, V = where, 

2B{6— 18a+6u;Q+6cja7— 187 

— 3(3+37+a) 2 (l+7) ' 

Also, the deceleration parameter reduces to, q = — 1§ = 37+ 2 " +1 < (using equation (21)) 
Hence, the present Universe is in a state of expansion with acceleration. 
Also, we get to = -Mill) _ 3±a and Q = _ 3( 1 + 2 t) 2 . 

Also, 7>— l^a<2 and w > — |. For the present Universe (i.e., taking 7 = 0) and the ACDM model, 

, , _ 3+a 

Case II: 

Now we choose u((j>) to be dependent on 0. Again we consider the power law forms, (12) and (17). Solving 
the equations in a similar manner, we get 



a8 + 2a + 3 - 8 1 + 7 „, , 3a ^+i)-2 3<*( 7 +i)+<3-2 2 A; , 2 < 1 -°) 

^Poao- 3(1+7) 0o 3 ^ ^ + io^^^ ( 22 ) 

and 



U(</>) = (2a + /?)(3a + 2/3 - 1)^,7 - (1 - 7 )poao~ 3(1+7) 0o^ ±1 V 2fi ^ ±il + -^U^^o^ (23) 
Substituting these values in equation (8), we get 

either (3 = -2 or 8 = -2a (24) 
Therefore for cosmic acceleration q < =>■ a > 1 and /3 < — 2. 
Therefore for the present era, 



G 




Fig. 3 and 4 shows the variation of u> against <f> for respectively closed and open models of the Universe. We take 
different values of a = 1, 1.5, 2, 2.5, 3. In both the figures we have considered the present dust filled epoch, i.e., 7 = 

and (3 = — 2, normalizing the parameters as a = po = 4>o = 1- 



u = ~- \ ^ and V = 2{a - l)(3a - 5)cj) a t - poa - 3 if = -2 

2 40o 



^, = _3 _ PoaoJ^ and = _0 nan -3 t i if p < _ 2 (25) 
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Also for vacuum dominated era, 



u = ~ and V = 2(a- l)(3a- 5)i - 2p for /3 = -2 
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w = -- and V = -2p for /3 < -2 (26) 



IV. MODEL USING GENERALIZED CHAPLYGIN GAS 



Here we consider the Universe to be filed with Generalized Chaplygin Gas with EOS 



Here the conservation equation (9) yields the solution for p as, 

C 



a 3(l+n) 

where C is an integration constant. 

A. Solution without potential: V{4>) = 

Case I: 



Fig.5 

V = V(<f>),P = -2,k = 



Fig.6 

barotropic fluid: (dust) 
V = V(<f>),P = -2,k = l 



Fig. 7 

V = V(<f>),P = -2,k = -l 



Fig. 5, 6 and 7 shows the variation of V against <j) f° r respectively flat, closed and open models of the Universe. We 
have considered different values of a = 1,1.5, 2, 2.5, 3 and f) = — 2. In all the three the figures we have considered the 
present dust filled epoch, i.e., 7 = 0, normalizing the parameters as do = po = (f>o = 1. 
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Fig.8 Fig.9 
barotropic fluid: (dust) 
V = V{4>),/3= -2a, 7 = w = w(0),ft = l,/3 = -2a, 7 = 



Fig. 8 shows the variation of V against the variation of <j> for all the models of the Universe, whereas, fig. 9 shows the 
variation of lu for only the closed model of the Universe. We have considered different values of a = 1, 1.5, 2, 2.5, 3 and 
the present dust filled epoch, i.e., 7 = 0, normalizing the parameters as ao = po = <t>o = 1, also as the calculation shows, 
for this (5 = —2a. For figure 8 the results for different values of a coincides with each other in each model of the 
Universe. 



First we choose uj(4>) = co ^constant. 
We consider the power law form 

= faa a (29) 

Equations (6), (7), (8) give, 

a 2 6 

{2uoa — 6)a + (uua 2 + Auoa — 6) — = -k (30) 



which yields the solution, 
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Fig.10 Fig.ll Fig.12 

Generalized Chaplygin Gas: 
lu = Lo(4>),f3 = -2,k = w = w(0),/3 = -2,fc = l w = = -2, fe = -1 



Fig. 10, 11 and 12 shows the variation of lu against <p for respectively flat, closed and open models of the Universe. We 
have considered different values of a = 1, 1.5, 2, 2.5, 3 and /3 = — 2. In all the three the figures we have considered n = 1, 

normalizing the parameters as ao = po = <f>o = B = C = 1. 



where P = — + 4 "° 6 and K is an integration constant. 

First we consider P > 0. Multiplying both sides of equation (31) by a p after squaring it, we get K n = 0, 
therefore giving, a = {uJ ^ 3)p t. 

Hence for flat Universe, we get, a =constant. 

For open model, we must have tua < 3 and a — \J ^ 3 ^ a -jp t, whereas, for closed model, uua > 3 and 
t. In all cases q = 0, i.e., we get uniform expansion. 



(W-3)P 

If P = 0, ad = ^k, i.e., a 2 = A In a + K . 

If k = 0, a = a/Ko^ + Co, (Co is an integration constant) causing g = 0, i.e., uniform expansion again. 
Case II: 

Now we consider u> = w (</>), i.e., u) dependent on </>. 

Also the power law forms considered will be (12) and (17). Solving the equations we get, 

, n 3a(l + n)-2 - 3c ( 1 + n) - + 2 2(l-c) 

W M = ^2 r 3a a + „) -,„, 1+ „, 1T fe + nmtEil ^ 



ft 2 



3a(l+n) -3a(l + n) 

5 + Cao-^+^o^ 3- 3 



Also substituting these values in the given equations, we get, either n = —1 or B = and also k = 0. If n = — 1, 
we get back barotropic fluid, and if B = 0, we get dust filled Universe. In both the cases the Generalized 
Chaplygin gas does not seem to have any additional effect on the cosmic acceleration. 



B. Solution with potential: V = V{<t>) 

Case I: 

Let us choose uj{4>) = lu =constant. 
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Fig.13 Fig.14 Fig.15 

Generalized Chap ly gin Gas: 
V = V(4>),(3 = -2,k = V = V{(j>),P = -2,k = 1 V = V(<f>),(3 = -2,k = -1 



Fig. 13, 14 and 15 shows the variation of V against <j) for respectively flat, closed and open models of the Universe. We 
have considered different values ofa = l,1.5,2,2.5,3 and (3 = — 2. In all the three the figures we have considered n = 1, 

normalizing the parameters as a = po = 4>o = B = C = 1. 



We again consider the power law forms (12) and (17). We get the solution for V((f>) to be 

p-2 -2B - Cao' 3(1+ " ) ii)n 3 °'g + °V 3 °' jl+ ''' 4fc 0-20 2c 
F(0) = (2a + P)(3a + 2/3 - 1)^^ + 2 ?° ^ + ^-^tt ( 33 ) 



3a(l + T») -3o(l + n) 

B + Ca - 3 ( 1 +")</)o^ 3— 3 



(i+n) a z 



Substituting these values in the other equations we get that n = — 1, i.e., the equation of state of Generalized 
Chaplygin Gas takes the form of that of barotropic fluid. Also we get, a = 1, which implies q = 0, i.e., uniform 
expansion of the Universe. 

Case II: 

Now we choose uj(</)) to be dependent on </>. 

Again we consider the power law forms, (12) and (17). Solving the equations we get the solutions for 
Brans-Dicke parameter and self-interacting potential as same as equations (32) and (33) respectively. 

Substituting these values in equation (8), we get 

either ft = -2 or 0= -2a (34) 
Therefore for cosmic acceleration q < => a > 1 and /3 < — 2. 
Therefore for the dust dominated era, 

3q-4 - w„ _ . 3o 



3 poop 3 (/>— 2(o-1)(3q-5) 2 _ 3 0V 
3^^ — an " ^ = — Poao — 55- 1/ P = — 2 



W =-0~ 3c-2 



3 _ jOpflp 3 <; 
2 4a 2 o ~ 



*-2 

-3 



and V" = -p ao if /? < -2 (35) 
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Also for vacuum dominated era, 
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Fig. 16, 17 and 18 shows the variation of V against <f> for respectively fiat, closed and open models of the Universe. We 
have considered different values of a = 1, 1.5, 2, 2. 5, 3 and (3 = — 2a. In all the three the figures we have considered 
n = 1, normalizing the parameters as ao = po = 4>o = B = C = 1. 
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a=2 
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Fig.19 Fig.20 
Generalized Chaplygin Gas: 
uj = oj(<j>),k = 0,j3= -2a cu = uj((j)), k = 1, (3 = -2a 



Fig. 19 and 20 shows the variation of u> against <j> for different values of a = 1, 1.5, 2, 2.5, 3. Here we have considered 
n — 0, normalizing the parameters as ao = po = 4>o = b = C — 1, also as the calculation shows, for this (3 = —2a. 



3 S 1 
u = -- and V = 2(a-l)(3a-5)^- -2\p vac ]„_ 2 for (3 = -2 
^ 00 
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w = -- and ^ = -2 [Pt,ac]2a+/3=o for < -2 (36) 



V. CONCLUSION 



We are considering Friedman-Robertson- Walker model in Brans-Dicke Theory with and without potential 
(V). Also we have considered the Brans-Dicke parameter (w) to be constant and variable. We take barotropic 
fluid and Generalized Chaplygin Gas as the concerned fluid. 

Using barotropic equation of state, we get, (i) for V = and ui ^constant, u> < and q < for some 
values of a, giving rise to cosmic acceleration, (ii) for V — and u> — uj{4>), we obtain cosmic acceleration 
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depending on some values of a and 0. In this case we get acceleration for closed model also at the radiation 
phase. We can show the variation of lo((/>) against the variation of <j> here [figure 1,2]. Figure 1 shows that 
as the value of a increases lo decreases steadily against the variation of <f>. For a > 1, we have accelerated 
expansion. The figure shows that the greatest value of lo can be — | and it decreases further as <p increases, 
(Hi) for V = V((f>) and lo ^constant, we get acceleration in the flat model irrespective of the values of a, (iv) 
for V = V(4>) and lo = u (</)), cosmic acceleration is obtained for (3 < 2. Here we can represent the variation 
of lo and V against the variation of <j> for (3 = —2 and (3 = —2a. For (3 = — 2, the variation of lo against <j> 
is same as figure 1 and that for closed and open models are given in figure 3 and 4. Here wc can see that 
for open model lo starting at — § decreases further, whereas for closed model lo starting at — | increases to 
e positive for a = 2. Figures 5, 6 and 7 sow that variation of V against the variation of <fr for (3 — —2 in 
respectively flat, closed and open models of the Universe. Here we can see that only for the closed model 
the potential increases positively, in the other two cases the potential becomes negative after a certain point. 
Figure 8 shows the variation of V against <j> for (3 = —2a. Again positive potential energy is obtained for only 
the closed model. The variation of to is shown in figure 9 for k = 1 and we can see that lo increases starting at — | . 

Using Generalized Chaplygin Gas , we get, (i) for V — and lo ^constant, uniform expansion is obtained, 
(ii) for V = and lo = lo(4>), Generalized Chaplygin Gas does not seem to have any effect of itself, (in) for 
V = V(4>) and to =constant, we get q — giving uniform expansion, (iv) for V — V(cf>) and lo = lo(4>) cosmic 
acceleration is obtained for (3 < 2 as previously obtained for barotropic fluid. Figures 10, 11, and 12 show 
the variation of lo for (3 = — 2 in flat closed and open models and the natures of the graphs do not vary much 
from that for barotropic fluid. Figures 13, 14 and 15 show the variation of V for flat, closed and open models 
respectively. Here for open model we get a negative potential after a certain point, whereas for closed model 
we get a positive potential always. For spatially flat model a positive V is obtained for a = 1.5, 2. Figures 16, 
17 and 18 show the variation of V for the models of the Universe for (3 = —2a. Positive potential is obtained 
for closed model and flat model shows positive potential for a > 1. For open model we get negative V again. 
Figures 19 and 20 show the variation of lo for flat and closed models respectively (f3 = —2a). For flat model to 
starting at — | decreases further and for closed model it increases slowly from — |. 

We have used B-D Theory to solve the problem of cosmic acceleration. Here we use barotropic fluid and 
Generalized Chaplygin Gas. Although the problem of fitting the value of lo to the limits imposed by the solar 
system experiments could not be solved fully, for closed Universe and (3 — —2 and a > 1, to starting from — | 
increases and for large <j>, we get to > 500, for both barotropic fluid and Generalized Chaplygin Gas. Also for 
flat Universe filled with barotropic fluid taking lo ^constant and V = V(<fi), we get the Bertolami-Martins [11] 
solution, i.e, V — V(<p 2 ) and q — — \ for a = At 3 . But taking Generalized Chaplygin Gas, we get accelerated 
expansion only when both lo and V are functions of the scalar field (f>. For (3 — —2 we get cosmic acceleration in 
the closed model, whereas, (3 = —2a gives acceleration in both closed and flat models of the Universe, although 
for flat Universe lo varies from — | to —2 and for closed Universe to takes large values for large 4>. In the end 
we see that for all the cases accelerated expansion can be achieved for closed model of the Universe for large 
values of to. Also the present day acceleration of the Universe can also be explained successfully, although in 
this case lo cannot meet the solar system limits. 
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